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Abstract. We show that the position operator in a class of /-deformed oscillators 
has a fractal spectrum, homeomorphic to the Cantor set, via a unitary transformation 
to Harper’s model. The class corresponds to a choice of ergodic operators for the 
deformation function. Hofstadter’s butterfly is plotted by direct diagonalization of a 
position operator with an originally vanishing diagonal. This is equivalent to a one¬ 
dimensional hamiltonian without potential. 
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1. Introduction 

The idea of deforming the oscillator algebra to obtain new solvable systems - in particular 
those connected to quantum optics umi - is a vivid example of an old precept: there 
are two classes of problems, given by the unsolved ones and the harmonic oscillator; 
the idea is to move examples of the former class into the latter [5]. In this respect, the 
q- and /- deformations of the Heisenberg algebra have proved to be sufficiently 
rich, containing hamiltonians with both finite and infinite levels in the spectrum [8]. 
Moreover, the inherent nonlinearity of deformations suggests a strong relation to old, 
but relevant models of electrodynamics [MI]. In a different realm of physics, but 
following the same precept, a mapping of a Dirac oscillator to a lattice has been reached 
by means of recursion relations This indicates a direct connection between the 

number of quanta of an oscillator and the site number of a semi-infinite chain [T3] . 

Despite the vast literature related to this subject, little has been said about 
the configuration space emanating from algebraic deformations. In this paper we 
show the existence of an /-deformed Heisenberg algebra that turns the hamiltonian 
H = + X'^) into an ergodic operator, while its position X - equivalently P - 

acquires a spectrum given by the Cantor set. This surprising result makes manifest the 
non-trivial topological properties of an algebraically modified space. For instance, the 
spectral measure of such a space vanishes, despite the uncountability of the eigenvalues 
associated to the deformed position operator. We thus have a purely singular phase 
space, if we borrow the term from spectral theory jl4j . 
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We proceed in two steps. In section H] we show a direct relation between the 
position operator of a deformed Heisenberg algebra and a deformed tight-binding chain 
without on-site potentials. Then, in section |3] we map unitarily an off-diagonal operator 
representing X to Harper’s model. The map contains as a special case a deformed semi- 
inhnite tight-binding chain without potentials. The spectrum of the resulting system 
is found to be identical to Hofstadter’s butterfly [15] when plotted as a function of an 
external parameter. A few remarks on hnite size effects in numerical results are given. 
We conclude in section 01 

2. A connection between the position of /-oscillators and a tight-binding 
chain without diagonal 

The /—deformations of Heisenberg’s algebra suggested by Man’ko [3] allow the 
introduction of a very general deformation function F{N) of the usual number operator 
N = ^{p‘^ + — 1). In natural units, we have the well-known relations 


[x, p] = / [a, a)~\ = 1, aF{N) = F{N -|- l)a. 


( 1 ) 


a = —-|- ip), = —;=(x — ip). 

y/2 V2 


( 2 ) 


Meanwhile, the deformed operators A, are establishecU by means of F{N) (hermitian) 


as 


A = F{N)a, A^ = a^F\N) = a^F{N), (3) 

with the following modihed commutations relations 

[A, = F{N)\N + 1) - F{N - 1)^N. (4) 

A deformed position X and momentum P can be constructed in the usual manner, but 
employing now A, Ah 


X = F{A + Af), P 

and the new commutation relation is 


1 

iy/2 


(A-At), 


( 5 ) 


[X, P] = i [F{Nf{N -Fl) - F{N - 1)2a] . (6) 

Evidently, these deformations have the effect of changing the spectrum of a hamiltonian 
conventionally dehned as a quadratic form in phase space. We have 


H = = i(AAt + At A) = i [F(A)2(Ar + i) + f{N - l)^^] , (7) 

f For convenience we use (H, which is slightly different from the standard notation. See the third 
relation in O- 
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showing that the spectrum cr^H) is not necessarily equispaced. Since the last relation 
expresses hT as a function of N, u{H) can be obtained by acting directly on the usual 
harmonic oscillator basis {|n)}n=o,i, 2 ,..., resulting in the following expression 

En = ]^[F{nf{n+l) + F{n-lfn\, n = 0,1, 2, 3,... (8) 

On the other hand, the spectrum of the position operator X can be found by 
diagonalizing the matrix representation {n\X\m) between harmonic oscillator states: 


{n\X\m) 


1 

7 ! 


F{n)\/n+ I5n,m-i + F{n 


l)\/n6n , 771+1 


Incidentally, since phase space operators are unitarily related by 


( 9 ) 


P = exp 



X exp 





( 10 ) 


then P and X are isospectral. Now, if the eigenstates |.^) of X are expanded in the 
usual oscillator basis |0 = 4>n\n) , we are led to a recurrence relation for cfn after 

substituting in X|,^) = i.e. 


F{n)\/n + 1 , F{n — l)y/n 

-<Pn+i + - 


1—1 


( 11 ) 


V2 +2 

with the boundary condition 0_i = 0. This can indeed be interpreted in terms of 
a nearest-neighbour tight-binding model; if we identify F{n — l)^yn/2 with hopping 
amplitudes (also known as couplings) and f with the k-th energy of a hopping 
particle in a chain, we have the equivalent relation 


= e.i'l ( 12 ) 

and the boundary condition can be hctitiously imposed by setting Aq = 0; this completes 
the analogy with a semi-inhnite chain. We may employ translation operators T, and 
a site-number operator N to express the tight-binding equation in hamiltonian form. 


Ht,,\ct>^) = [A(A)T + T^X\N)] |+) = e^\ct>^), (13) 

OO 

i/) = E'^»i’i>' (14) 

n=0 

with the properties 

T|n) = |n-t-1), T^|n) = |n —1), T'^|0) = 0, N\n)=n\n). 


(15) 
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The atomic states \n) must now be interpreted as vectors whose space representations 
{y\n) are functions of a real position variable y; such functions are localized around site 
n = 0,1,2,.... It is convenient to mention here that in certain applications of tight- 
binding models, such as photonic crystals [16] or optical lattices [I71[T8], the coupling 
constants A can be engineered by means of evanescent waves coming out of a potential 
well or a resonating cavity located at a lattice site. 

2.1. Ergodic operators for deformed oscillators 

Ergodic hamiltonians [I9] are a key ingredient of fractal spectra, in particular for the 
‘almost’ Mathieu operator [20]. In a more general fashion, an ergodic hamiltonian may 
contain operators whose spectrum is dense and uncountable in a subset of M. Here we 
must specify our deformation function such that + A^) is ergodic. 

An operator O with such a property can be obtained by means of a quasiperiodic 
function of N. Let O be periodic of period r: 


d{N + T) = d{N) (16) 

but now, if u and r are non-commensurate, the operator 0{uN) does not have period 
T in the variable N: 


0{cj{N + t)) = 0{ojN + cor) 0 {ojN). (17) 

Moreover, there is no translation of oscillator quanta that leads to invariance, i.e. for 
any k,q E Z 


0{uj{N + k)) = OiujN + koj) ^ 0 {ujN -|- qr) = 0{uN). (18) 

In addition, if O is a hermitian bounded operator its spectrum must be dense, for acting 
on In) yields a set equivalent to the irrational numbers within [Omin,Omax]- A more 
specihc example can be given in terms of trigonometric functions: If our deformation 
function F is set as 


^ smlu,{N+l)] 

^ ’ x/aTT 

then the hamiltonian H = + A^) will become 

H = ^ (sin^ [oj{N + 1)] + sin^ [wA]) . (20) 

Now we take uj to be non commensurate with tt, such that En generates a dense set 
contained in an interval. One can verify that the limits of the interval depend on u. 
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3. A map between position operators and Harper’s hamiltonian: the 
realization of Hofstadter’s bntterfly 

The position operator X now has the form of an ergodic operator, since it contains ffT^ 
in its definition, and when u is non commensnrate with vr, the spectrnm of sm{uN) is 
dense in the interval [—1,1], Here we claim and prove that in fact X is isospectral to a 
tight-binding hamiltonian that prodnces Harper’s eqnation, namely 


-f^Harper = T 2 COS {uN - I/) . (21) 

This is in principle nnusnal, since -Hnarper contains no lattice deformations (constant 
conplings) and the effect of fractality in its spectrnm is due solely to the external field. 
On the other hand, X is a pure deformation, without diagonal terms. We shall see that 
there exists a unitary U such that UXU^ = ATTnarper with A an overall scaling factor. 
Let us start with 


^ sin[a;(X + l)] tsin[w(X + l)] 

j\ — . (Jj (Jj . 

V‘^{N + i) VWTT) 

sin [wX] ^ sin [wX] ^ 

— a / “t” / a 

cos [cjX — 7r/2] ^ cos [caX — 7r/2] .|. 

— a . , a . 


( 22 ) 


From here, we see that we may work with more generality using the following rescaled 
operator and the introduction of an extra parameter v 


X„ = a 


2 cos [cjX — z/] ^ 2 cos [caX — \ 


y/N 


y/N 


exp [i{u;N — z/)] ^ exp [—i{uN — z/)] .j. ^ exp [—i{uN — z/)] ^ exp [i{uN — z/)] ^ 


= a 
= Ho+ V, 
where 


y/N 


y/N 


a' -|- a- 


y/N 


+ 


y/N 


(23) 


Ho = a 


exp [i{ojN — v)] ^ exp [—i{ujN — v)] ^ 


Vn 


y/N 


exp [—iipjN — i/)] exp [i{uN — z/)] .|. 
\/X y/N " 


V = a 


a'. 


(24) 


In what follows we show that Ho corresponds to the addition of two translation operators 
and H to a local potential of a transformed variable, using the appropriate basis. In 
order to transform X^ unitarily, we work now at the level of states using the definition 
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with J\f an appropriate normalization constant. It is a matter of simple algebra to verify 
the following equations: 






^ exp - 01 exp[.(^-0] ^q ^ 


(26) 


(27) 


Here one may proceed by direct substitution of fl25p in the second relation of fl26p . while 
the hrst relation follows by recognizing that 


-/^XPH^-Ol^tjm. 


(28) 


In fact, the expansion coefficients in fl25l) can be derived by solving the set of recurrences 
emerging from fl26l 1271) as requirements. Now that we have established these properties, 
the following identifications are in order 


Ho = Tk + tI, V = 2 cos [uK — v] (29) 

i.e. Tk,Tl are translation operators in k and H is a function of the number operator K 
such that K\k) = k\k). We also write a) = 2a;, k = 2v io simplify the notation. Harper’s 
equation now emerges transparently when Xy is applied to its eigenvector |^) written 
as a linear superposition |^) = YlT=-oo'^k\k)'^ recurrence is 


fjk+i + ^Jk-iF 2, cos {uk - ■u)tlJk = (30) 

Moreover, since /c G Z, this equation corresponds to an inhnite chain. 

It is left to verify that the transformation \n) i—)■ |/c) is unitary. This is most easily 
done by applying the transpose conjugate of fl25l) to \k) as follows: 


^ exp {— i [a; (|n(n + 1) + k{k — 1 — 2n)) + z/n] } \k) 


k=—oo 
- oo 


jjr ^ exp [—i [a; + 1) + k{k — 1 — 2n)) + z/n] } 


k=—oo 


X jjr exp [a; (im(m + 1) + k{k — 1 — 2m)) + z/m] } |m) 


m=0 


exp [i [|a;(m(m + 1) — n{n + 1)) + z/(m — (f^)] } ^ 


^2iuk{n—m) 


m=0 


\m) 


k=—oo 


= |n) 


(31) 
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where in the last line we recognize that the summation in braces is a properly normalized 
delta, so m = n. 

It is interesting to note that the map described above does not work when (1^ 
contains the more general potential 2Acos{u}N — u) with A 7 ^ 1. This seems to be 
connected with the known result na that only A = 1 produces a purely singular 
spectrum. For general A, the deformation function becomes complex 


F{N) 


exp {i [cc;(A^ + 1) - + Aexp{-i [a;(A^ + 1) - i^]} 

\/FTT 

+ A2 + 2A cos [2u:{N + 1) - 2z/] 

x/FTI 


(32) 


and its spectrum is modified. When A S> 1 or A <C 1, the trigonometric part becomes 
a mere perturbation. 


3.1. The spectrum of X,y 

As pointed out in very early works [T51I2I] on this subject, the spectrum of this system 
is equivalent to Cantor’s set when u is non commensurate with vr, i.e. it is purely 
singular. This is independent of the value of z/. On the other hand, if u and vr are 
commensurate we have qoo = pvr, and the system becomes periodic. The period is q if 
p, q are relatively primes; this implies that a structure of q bands should emerge. As 
p/q approaches irrationality, the number of such bands increases as a power of a given 
initial number (say, two bands for a precursor of Cantor’s set) and their length decreases 
in a similar amount. This analysis is also valid for our off-diagonal Xy as expressed in 
fl2^ . since it only requires the ergodic behavior of sin(a;A^ -|- v). Usually, the spectrum 
is depicted by showing a sequence of shrinking intervals as a function of ca, giving rise to 
the famous butterfly. In our problem, we may diagonalize directly the operator X^ by 
a suitable numerical routine and plot the spectrum. Since the matrix representations 
that we use are of a hnite size, an additional artefact appears in the pattern: The 
gaps between bands for rational u contain eigenvalues due to edge states. Such bands 
proliferate in a self-similar pattern and it is to be expected that edge states appear in 
various regions of the spectral span. To verify that indeed such states are due to hnite 
size effects, we proceed to diagonalize X^, for increasing sizes rimax = 200, 210,... and we 
vary u = 7r/200, 27r/200, 37r/200,... for each size. We show the results in £g. [H In £g. 
[ 2 ] we repeat the calculation, but with periodic boundary conditions; the results reveal a 
signihcant displacement of border levels towards the allowed bands. 


4. Conclusion 


A fractal position spectrum has been found. The example that we have studied has been 
shown to be equivalent to the reduced problem of an electron in a square lattice under the 
inhuence of a magnetic held. We also computed the spectrum numerically to verify the 
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Figure 1. Numerical spectrum of two different off-diagonal ergodic operators X^. 
Right panel: rimax = 200. Left panel: nmax = 210. The red arrows indicate a 
variation of edge states according to the size of the arrays. In both cases, the abscissa 
is partitioned in 200 fractions of tt. 



Figure 2. Numerical spectrum of the ergodic operator X^, rimax = 200 and periodic 
boundary conditions. The red arrow shows the location of edge states, at their new 
position. 


property of isospectrality, recognizing that very small discrepancies arise dne to finite 
size. Undonbtedly, experiments can be carried ont in the context of artihcial solids, 
snch as those produced by optical lattices, photonic crystals or microwave resonators. 
In particular, a nearest-neighbour tight-binding chain of identical sites and adjustable 
distances has been used succesfully to produce other models, and this choice is likely to 
give immediate results in the present case. It is also important to note that a realization 
of ergodic hamiltonians using a hxed photon mode in a cavity is not a complete diversion 
from quantum optical applications: it is a non-linear example of oscillation that can be 
studied and perhaps partially tested. Due to the extreme dependence of its effective 
frequency Q{n) = E^/{n+1/2), we may propose D to vary randomly with the occupation 
number. 
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